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ABSTRACT 

A  cloud- ocean  planetary  boundary  layer  (OPBL)  feedback  mechanism  is  presented  and  tested  in  this  paper. 
Water  vapor,  evaporated  from  the  ocean  surface  or  transported  by  the  large-scale  air  flow,  often  forms  convective 
clouds  under  a  conditionally  unstable  lapse  rate.  The  variable  cloud  cover  and  rainfall  may  have  positive  and  negative 
feedback  with  the  ocean  mixed  layer  temperature  and  salinity  structure.  The  coupling  of  the  simplified  Kuo's  (1965) 
cumulus  cloud  model  to  the  Kraus— Turner's  (1967)  ocean  mixed  layer  model  shows  the  existence  of  this  feedback 
mechanism.  The  theory  also  predicts  the  generation  of  low  frequency  oscillation  in  the  atmosphere  and  oceans. 

I.  INTRODUCTION 

The  feedback  mechanism  between  clouds  and  ocean  mixed  layer  can  be  illustrated  as  fol¬ 
lows.  First,  clouds  reduce  the  incoming  solar  radiation  at  the  ocean  surface  by  scattering  and 
absorption,  which  cools  (relatively)  the  ocean  surface  layer  by  increasing  mixed  layer 
entrainment.  The  cooling  of  the  ocean  mixed  layer  lowers  the  evaporation  rate,  which  will 
diminish  the  clouds.  This  is  a  negative  feedback  mechanism.  Second,  precipitation  dilutes  the 
surface  salinity,  stabilizing  the  upper  ocean  and  reducing  mixed  layer  deepening.  The  mixed 
layer  may  even  be  caused  to  shallow  if  the  downward  surface  buoyancy  flux  is  sufficiently  en¬ 
hanced  by  the  precipitation.  The  reduction  in  mixed  layer  depth  will  increase  the  sea  surface 
temperature  (SST)  by  concentrating  the  net  radiation  plus  heat  fluxed  downward  across  the 
sea  surface  into  a  thinner  layer.  The  increase  of  SST  augments  the  surface  evaporation,  which 
in  turn  produces  more  clouds.  This  is  a  positive  feedback  mechanism.  Fig.  1  shows  the  main 
physical  processes  (heat,  mass,  and  momentum  fluxes)  at  the  two  adjacent  boundary  layers: 
the  OPBL  and  the  marine  atmospheric  boundary  layer  (MABL).  Fig.  2  illustrates  this 
positive  /  negative  feedback  mechanism. 

Since  clouds  have  significant  effects  on  the  large-scale  atmospheric  circulation  through 
the  transfer  of  heat,  moisture,  and  momentum,  and  on  the  ocean  mixed  layer  through  the 
attenuation  of  solar  radiation  at  the  ocean  surface,  and  since  the  SST  is  an  important  factor 
for  the  development  of  clouds,  the  feedback  mechanism  mentioned  above  has  a  potentially 
significant  impact  on  the  air— sea  interaction,  weather  and  ocean  prediction. 

The  coupled  model  is  one-dimensional.  We  are  aware  of  the  importance  of  the  horizon¬ 
tal  advection  and  the  limitation  of  one-dimensional  models.  However,  the  intent  of  this  work 
is  to  develop  a  formalism  to  examine  this  thermodynamic  feedback  between  the  two  fluids. 
Because  we  wish  to  concentrate  on  the  thermodynamic  interaction,  horizontal  advection  is  ig¬ 
nored  initially. 
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Fig.  1.  Main  physical  processes  at  the  two  adjacent  Fig.  2,  Feedback  mechanism  between  OI’Bl  and 

boundary  layers.  clouds. 

II  THE  OCEANIC  MIXED  LAYER 


I .  Basic  Equations 

Much  of  the  one-dimensional  theory  for  the  OPBL  or  mixed  layer  is  dependent  upon  the 
validity  of  two  crucial  hypotheses.  The  first  of  these  is  that  vertical  mixing  within  the  turbu¬ 
lent  boundary  layer  and  entrainment  mixing  at  its  base  occur  in  response  to  the  local  atmos¬ 
pheric  forcing-  the  surface  wind  stress  and  the  buoyancy  flux  at  the  sea  surface.  The  second 
hypothesis  is  that  the  mechanical  energy  budget  is  the  key  to  the  understanding  and  predic¬ 
tion  of  mixed  layer  dynamics  (Garwood,  1979). 

The  buoyancy  flux  is  attributable  to  heat  flux,  evaporation,  and  precipitation,  and  the 
shear  production  of  turbulence  is  attributable  to  surface  wind  stress.  The  mixed  layer  temper¬ 
ature  T„  salinity  S„  and  depth  h  are  predicted  by  a  simplified  form  of  the  Garwood 
(1977)  mixed  layer  model.  It  is  essentially  a  “calibrated’  Kraus  and  Turner  (1967)  model 

which  is  modified  to  include  salinity  and  advection: 

cT  F  j 

h~jr  =  ~wATs  -r-y,) — >  (D 

Cl  pc 

r  w  p 

es 

h-£-=-W'(S!-S_h)  +  (E-Pr)S!,  (2) 

where  h  is  the  mixed  layer  depth,  c is  the  water  specific  heat  under  constant  pressure, 
pn  is  the  characteristic  water  density,  E  is  the  surface  evaporation  (m  /  s),  Pr  is  the  precip¬ 
itation  rate  (m  /  s),  and  w,  is  the  entrainment  velocity  parameterized  as  (Chu  and  Garwood, 
1988): 

(C,u\  -  C2B0h ) 

"  gf>HTl-T,t)-^Ss  -5.,)]  ’ 
where  u ,  is  the  water  surface  friction  velocity,  which  is  computed  by 


(3) 
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where  CD  is  the  drag  coefficient,  taken  as  0.001,  |  V*"’  |  is  the  wind  speed  at  a  10  m  height, 
pa  is  air  density  near  ocean  surface,  and  p„.  is  the  sea  water  density.  For  a  wind  speed  of 
10  m  /  s,  u ,  —  1 , 1  cm  /  s.  For  the  case  of  an  ocean  surface  without  ice  (low  and  middle  lati¬ 
tudes)  the  surface  buoyancy  flux,  flf„  has  two  components: 

*g  F  I 

50= - (5) 

Here  a  is  the  water  thermal  expansion  coefficient,  and  fi  is  the  salinity  contraction 
coefficient. 

The  effects  of  clouds  on  the  buoyancy  flux  at  the  ocean  surface  are  two-fold:  (1)  de¬ 
creasing  through  the  increase  in  the  net  heat  loss  at  the  ocean  surface,  Fr  |  ,  by  re¬ 

ducing  the  incoming  solar  radiation,  and  (2)  increasing  Ba  due  to  precipitation. 

The  surface  heat  flux,  Fr  I  (upward  positive),  is  computed  by 

'  'l-o 

=  <*) 
Here  Rs  is  the  incoming  solar  radiation  absorbed  by  the  ocean  surface,  Rb  is  the  net  ener¬ 
gy  loss  from  the  ocean  surface  through  longwave  radiation,  L  is  the  latent  heat  of 
vaporization  of  water,  H,  is  the  sensible  heat  flux  to  the  air.  The  variables  T_h  and  S.k 
are  the  temperature  and  salinity  of  the  water  immediately  below  the  mixed  layer  that  is  to  be 
entrained  into  the  mixed  layer.  The  standard  bulk  formulae  are  used  to  calculate  the  surface 
evaporation: 

E  =  PaCy:)\lqs(T>)-q,]/pit  (7a) 

and  the  sensible  heat  flux  from  the  ocean  surface: 


(7b) 

where  i/,(T)  is  the  saturated  mixing  ratio,  and  Ta0  is  the  air  temperature  near  the  ocean  sur¬ 
face. 

The  time  rate  of  change  in  mixed  layer  depth  is  computed  by 

^  =  w  —  w  ,  (8) 

ol  '  -* 

where  w-h  is  the  upwelling  velocity  at  the  bottom  of  the  mixed  layer.  In  order  to  slove  these 
equations,  we  should  compute  each  term  in  the  surface  heat  and  salinity  fluxes. 

2,  Cloud  Effects  on  the  Net  Radiation  at  the  Ocean  Surface 

Clouds  reduce  the  solar  radiation  upon  the  ocean  surface  by  scattering  and  absorption, 
which  is  computed  by  Budyko's(1978)  formula 


Here  Rb0  (340  W  ttT’)  is  the  solar  radiation  absorbed  by  the  ocean  surface  layer  under  a 
clear  sky.  The  parameters  xM  and  aa  represent  albedos  of  the  earth-atmosphere  system 
with  complete  cloud  cover  and  a  cloudless  sky  respectively,  and  have  the  following  values: 


x  =0.46,  a.  =0.2. 

SR  ’  iC 

The  ocean  surface  emits  longwave  radiation  to  the  atmosphere  and  to  space.  However, 
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clouds,  as  well  as  dry  air,  partially  absorb  the  radiation  and  re-emit  longwave  radiation  back 
to  the  ocean  surface.  Thus  the  net  upward  energy  loss  by  longwave  radiation  at  the  ocean  sur¬ 
face,  R„,  is  corrected  for  the  downward  radiation  by  the  clouds  and  the  air.  From  longwave 
radiation  data,  Budyko  (1978)  derived  a  semi-empirical  formula: 

Rt  =a  +  bT-(at  +h,)n.  (10) 

The  dimensional  coefficients  a,  b,  a,,  and  bx  are 

a  =  ~  377.6Wm  \  *  =  2.2  Wm  2K  ', 


a,  =  —  389 .8Wm 


b ,  =1.6  Wm 


III.  TIME  RATE  OF  CHANGE  OF  THE  CLOUD  COVER 

The  time  rate  of  change  of  cloud  cover  is  assumed  to  be  proportional  to  the  moisture 
supply  divided  by  the  amount  of  water  vapor  necessary  to  produce  the  model  cloud.  The  main 
processes  causing  the  cloud  dissolution  are  precipitation  and  mixing  with  the  environmental 
air.  The  cloud  evaporation  due  to  mixing  with  ambient  air  is  a  complicated  problem,  and  neg¬ 
lected  for  the  sake  of  simplicity  in  this  paper.  Thus  the  equation  for  cloud  cover  is 

on  (M,0+£-P) 

dt  W  '  K  J 

where  W  is  the  total  amount  of  water  vapor  needed  to  create  the  cloud  over  a  unit  area.  The 
large-scale  horizontal  moisture  convergence  in  the  column  of  atmosphere  per  unit  area  is  de¬ 
noted  From  mean  distributions  of  temperature  and  mixing  ratio  in  the  environmental 

air  and  inside  a  deep  cumulus  cloud  (Kuo,  1965),  we  estimate  that  W~  5  cm.  Eqs. 
(l),(2).(3),(8),  and  (1 1)  are  the  basic  equations  for  the  air-ocean  coupled  system. 

IV.  RELATIONSHIP  BETWEEN  PRECIPITATION  RATE  AND  CLOUD  COVER 

By  linear  regression  with  use  of  hourly  rain  amounts  and  satellite  IR  brightness  measured 
during  Phase  1,11,  and  III  of  GATE,  Albright  et  al.  (1985)  presented  a  linear  relationship  be¬ 
tween  average  precipitation  rate  P,  in  boxes  1.5°  (168  km)  on  a  side  and  cloud  cover  n 
of  the  boxes  by  clouds  with  tops  colder  than  —36)7 .  The  relationship  is 

Pr  (ms  ')  =  (0.472  +  8.333n)  x  1CT7.  (12) 

This  result  coincides  with  Arkin  (1979)'s  earlier  analysis  during  GATE  over  the  B— scale  ar- 


V.  MEAN  STATE  AND  PERTURBATIONS 


Mean  state  of  the  coupled  cloud-OPBL  system  (Ts,  Sr,  A,  £,  Pr,  )  is  obtained  by 
the  zero  time  rate  of  change  in  the  prognostic  equations  (1),  (2),  (8),  and  (11).  The  mean 
entrainment  velocity  w  can  be  obtained  from  (3) 


C7,»,  -C,Bah 


Ab=g(aAT-fiAS) 

is  the  reduced  gravity,  and  AT  and  AS  are  the  mean  temperature  and  salinity  jumps  at  the 
bottom  of  the  mixed  layer  _  _ 


AT=  T  ~T 


AS^S-S 


When  the  conpled  system  is  perturbed  from  its  mean  state,  the  thermodynamic  feedback 
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mechanism  between  the  cumulus  clouds  and  the  oceanic  mixed  layer  makes  the  perturbation 
either  grow  (positive  feedback)  or  dampen  (negative  feedback).  The  principal  purpose  here  is 
to  study  thermodynamic  feedback  mechanisms  between  clouds  and  the  oceanic  mixed  layer. 
Hence  the  energy  exchange  at  the  air-ocean  interface  is  a  primary  focal  point.  Therefore,  we 
shall  neglect  initially  the  perturbations  of  those  variables  which  are  not  directly  related  to  the 
exchange  at  the  air-ocean  interface.  Wyrtki  (1981)  estimated  that  the  mean  upwelling  in  the 
equatorial  Pacific  is  around  I  m  /  day.  Climatic  data  also  indicate  that  the  mean  evaporation 
and  precipitation  rates  are  on  the  order  of  1  m/  yr  (3.17*  Ifl'Vts-1).  The  direct  contribu¬ 
tions  of  evaporation  and  precipitation  to  the  change  of  mixed  layer  depth  are  neglected,  com¬ 
pared  to  the  contributions  due  to  upwelling  and  entrainment. 

The  perturbations  satisfy  the  following  equations: 


dr 

5 

~w 


=  -w  .r  -w'  AT-- 

-Ay  e 


F' 


p  c 

•  w  p 


_  ds'  -  - 

b  lf  =  -w^S'i-w'eAS  +  {E-Pr)S’s+S,(E'-P'r), 

dt  W( 
oh’ 

-J—  -  W  , 

dt  • 


— 3 

U  . 


w'A b=  -gv.Svr'-ps'j-c—h' 

h 


+  C, 


«gF, 


p  c 

r »  p 


ir*-  +fess(E'-p'r)+fe(E~pr 


(14) 

(15) 

(16) 
(17) 


(18) 


Eliminating  h‘ ,  w'e,  E',  P“r.  from  (14)— (18)  and  neglecting  the  small  terms,  we  get  the  fol¬ 
lowing  fourth  order  equation  for  each  variable: 

(19) 


,d"  53  .  s2  .  e 


where 


+av  +a’Tt+a<W'  n')  =  0’ 


<j>r=ps',~<xr„ 

_  i  c 


a^t+RH.+W 


_  m*. ,  1  , . 

a ,  =  -  — (C,  +  -)*  + 


(1  +  C,)U- 


2  C, 

C  ,mX 


+  (- 


'->f 


Rh.  W  ’ 


Riz. 


+  -2-0- 
TH-4T 


C,  mx  _T  ,, 


(20) 


Here  three  constants  t  ,  t  ,  r  are  time-scales  of  internal  motion,  ocean  turbulence,  and 
surface  evaporation  defined  by  _  _  _ 

_  h  _  h  _ h  y-j  i , 

t.  =— -  ,  T  ,  t  =-  ,  (21) 
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where 

i «  = 


L  ?E 

„r,1>  cT 


K 


dq,{Ts )  nc  ,  , 

=  0.5  x  10  ms  . 


P„c. 


dt 


The  constants  y  and  m  are  defined  by 

dP  I 

y=  !/f 


<»»> 

m  s-2—  =  8.5  . 
<xL 


Using  (1 1)  we  have 

Ri  is  the  Richardson  number 


y  =  8.333  x  10“  ms 


(22) 


(23a) 


(23b) 


Ris- 


Abr 


(23c) 


In  this  paper  Ri  is  taken  as  104.  The  two  model  parameters  A,  and  k  are  defined  by 

gW  On  ' 

which  is  the  dependence  of  the  surface  buoyancy  flux  change  on  cloud  cover,  and 


(24a) 


_  xgAT 
Ab 


(24b) 


which  is  the  fraction  of  the  mixed  layer  base  density  jump  induced  by  temperature  to  the  total 
density  jump. 

VI.  INSTABILITY  AND  OSCILLATION  CRITERIA 

The  general  solutions  of  (19)  have  the  following  forms: 

<p'(?)  =  Xy,.exp(<7  J),  «'(/)  =  £>.  export).  (25) 

where  djr  e,  (J=  1,  2,  3,  4)  are  the  integral  constants,  and  cr},  and  <r4  are  the 

eigenvalues  which  are  the  roots  of  the  fourth— older  algebraic  equation. 


4  3  2 

cr  +a,i t  +a2«r  +alu  +  aA=0. 


(26) 

For  the  purpose  of  a  preliminary  sensitivity  analysis  only  tc  and  A  are  allowed  to  vary,  de¬ 
pending  on  observations.  Other  parameters  are  given  constant  values  listed  in  Table  1  By  the 
definition  of  X  (24a)  and  the  estimations  of  y  (23b)  and  Budyko's  formulae  (8,  9)  we  can  es¬ 
timate  the  value  of  X  as 

2-0.45  x  kTV‘  (27) 

which  indicates  that  A  has  an  order  of  10“ V.  Therefore  in  this  paper,  A  varies  from 
-l<rV!  to  10“Vj.  It  is  reasonable  to  let  k  vary  between -10  to  10.  The  case  of  x  =l 
means  that  salinity  is  homogeneous  across  the  mixed  layer  base. 

The  instability  criteria  for  the  thermodynamically  coupled  air-ocean  system  are 

<0  decaying 

neutral,  (28) 

>  0  growing 

where  a  is  the  root  of  the  fourth-order  equation  (26).  The  oscillation  criteria  for  the  coup¬ 
led  system  are 

/m(?){=°  n°noSClUatory  (29) 

¥=■  0  oscillatory  . 

We  compute  all  roots  of  (26)  for  different  values  of  the  parameters  k  and  A,  and  ob¬ 
tain  four  roots  at  each  points  of  the  parameter  space  (ic.  A).  Two  roots  among  the  four  have 


Re(c 


_f  <0 

ffh  =0 

*  -  rs 
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negative  real  parts  throughout  the  whole  parameter  space,  representing  the  damping  modes, 
in  which  we  are  not  interested  here.  The  other  two  roots,  ert  and  er2,  have  positive  real 
parts  somewhere  in  the  parameter  space,  representing  the  existence  of  growing  modes  in  cer¬ 
tain  parameter  ranges. 

Fig.  3  shows  the  isolines  of  (a)  /?e(<r,)  and  (b)  | /w(<r x )|  in  the  (k,  A)  plane.  Fig.  4  in¬ 
dicates  the  isolines  of  (a)  Re(p2)  and  (b)  |/w(ct,)|  in  the  (k,  X)  plane.  Both  the  real  part 
(growth  rate)  and  the  absolute  value  of  the  imaginary  part  (periodicity)  of  the  roots 
a ,,  a 2  show  saddle-type  distributions.  Several  interesting  results  from  Figs.  3  and  4  are 
summarized  as  follows: 

(i)  A  necessary  condition  for  the  generation  of  growing  modes,  which  can  be  seen  from 
these  figures,  is: 

kA  >0.  (30) 

This  implies  that  if  the  surface  water  is  wanner  (cooler)  than  the  deep  water,  i.e., 
k>  0(k<  0),  the  growing  modes  are  excited.  In  this  case  that  the  surface  buoyancy  flux  in¬ 
creases  with  increasing  (decreasing)  cloud  cover. 

(ii)  If  the  damping  modes  are  ignored,  the  two  eigenvalues  <r1  and  flj  exhibit  the  simi¬ 
lar  properties.  The  growing  modes,  and  er2,  are  further  separated  into  oscillatory  and 
nonoscillatory  modes.  The  condition  for  this  separation  is  approximately  given  by 


kX 


{ 


<A  Oscillatory 
>  A  Nonoscillatory, 
where  A  is  a  positive  number  of  about  2  x  10'V . 

Combining  (30)  with  (31)  and  using  (24a),  the  necessary  condition  for  oscillatorily  grow¬ 
ing  modes  is 

K  O  —  »(/{,,  ~  R  ) 

A>WTn^S^~-^i  >° 


(31) 


(32) 


and  the  condition  for  the  nonoscillatorily  growing  modes  is 
k.  ft  _  xf/T  —  R  ) 

K  I  >A 


W  8n 1 


P  C 

r  It  p 


(33) 


The  properties  of  thermodynamic  instability  of  the  coupled  system  depend  largely  on  the  rela¬ 
tionship  between  precipitation  Pr  and  the  cloud  cover,  and  on  the  relationship  between  net 
radiation  at  the  ocean  surface  and  the  cloud  cover. 

(iii)  The  product  of  the  two  parameters  k  and  X  represents  the  relative  strength  of 
positive  to  negative  feedback.  When  kX  is  larger  than  the  critical  value  A,  the  positive 
feedback  greatly  exceeds  the  negative  feedback.  The  coupled  system  becomes  nonoscillatory 
and  growing.  When  kX  is  positive  and  smaller  than  the  critical  value  A,  the  positive 
feedback  slightly  exceeds  the  negative  feedback.  The  coupled  system  is  oscillatory  and  grow¬ 
ing.  When  kX  is  negative,  the  negative  feedback  exceeds  the  positive  feedback.  The  coupled 
system  becomes  damping.  Comparing  Fig.  4  with  Fig.3,  we  find  that  for  the  oscillatorily 
growing  modes  the  two  roots  have  the  same  growth  rate  and  frequency.  However,  for  the 
nonoscillatorily  growing  modes  the  growth  rate  relating  to  the  eigenvalue  tr2  is  much  larger 
than  that  relating  to  trt. 

(iv)  For  the  oscillatorily  growing  modes,  the  growth  rate  ar  has  the  order  of  0.5  * 
10'V1.  and  |or .  |  has  the  order  of  10 "V.  The  corresponding  period  for  the  oscillation  is 

T  =  rL\  =2yr.  (34) 
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(a)  Rc(«i)  (b)|tm(<r,)| 


Fig.  3.  Distributions  of  eigenvalue  al  in  the  ( k ,  >.)  plane  for  standard  case  (unit  in  10  ’s'1):  (a)  /feta,),  (b) 
|/'"(a1)|. 


(a)  Retaj)  (b)  |lm(<i2)l 


Fig.  -I  Distributions  of  eigenvalue  Oj  in  the  far,  l)  plane  for  standard  case  (unit  in  10  V1):  <a>  fle(a2),  (b) 
|/m(a,)|. 

VII.  SENSITIVITY  OF  |o.  |  TO  d  AND  * 

In  order  to  investigate  the  dependence  of  <r,  and  | erj  |  on  the  parameters  k  and  2,  we 
plot  rr,r  cr,,,  | o’ |.  jdj,  |  versus  k  at  five  different  values  of  k(k=  1,  5,  10,  -5,  —10),  and 
taking  the  same  values  for  the  parameters  as  in  Table  1 .  The  results  are  listed  below. 

(i)  Fig.  5  represents  the  case  for  k=  1,  i.e.,  the  case  with  no  salinity  jump  at  the  mixed 

layer  base.  The  identity  of  <rlr,  |trlf  j),  to  (<rir,  |a-a  j)  means  the  two  eigenvalues  are  com¬ 
plex  conjugates.  Fig.  5  shows  that  the  growth  rate  trr  is  a  linearly  increasing  function  of  k 
from  cr,  =  —0.5  X  10“V’  (damping  mode)  for  /.  =  —  10_6s~l  to  <Jr  =  0.5x  10_V  (growing 
mode)  for  2  =  0.5  x  10'V1.  The  frequency  |cr.  |  decreases  almost  linearly  with  k  from 
1(7,1=1.7x10  s  ’  (1.2  year  period)  for  2= -lO'V1  to  jo-,  |  =  0.86  x  10  1  (2.3  year 

period)  for  2  =  0.5  x  10~^s~'. 

(ii)  The  dependence  of  growth  rates  ulr,  er2r  and  frequencies  |<7|(  |,  |tr2(  |  on  2  for 
k  =  5  and  k=  10  is  depicted  in  Fig.  6.  The  value  of  <r2r  monotonically  increases  with  2, 
whereas  <jlr  has  a  single  peak  which  is  around  0.5  x  10_7s~’,  and  asymptotically  tends  to  ze¬ 
ro  as  2  ’■oc  The  necessary  condition  for  the  growing  oscillation  is 


-MH)  I  'o  | 


No.  1 
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0<H<Ac, 

A 

where  ,  is  a  critical  value  for 

r  K 


Fig.  5.  Dependence  of  {ay,  |er  N  |  and  (olr,  Fig.  6.  Dependence  of  (a j)  and  (cr^. 

|e 2i  ])  on  a  for  the  ease  without  salinity  jump  jor^  |)  on  J  for(a)  k=5,  (b)K=10. 

at  the  milted  layer  base  (x  -  1 ). 


Table.  1.  The  Values  of  the  Model  Parameters 


VIII.  SENSITIVITY  OF  a,.  |o]  TO  w_c 

To  understand  the  dependence  of  crr  and  |<x  |  on  the  mean  vertical  velocity  at  the  mixed 
layer  base,  iv  _  t ,  we  compute  the  roots  of  (26)  by  using  the  same  values  for  the  parameters 
as  previous  two  sections  except  for  w  _  t ,  which  is  changed  to  two  alternate  values:  w  h 
=- 0  2m /day  (weak  upwelling case)  and  w_li  =  2m  /  day  (strong  upwelling  case). 

The  main  results  are  that  the  growth  rate  and  the  frequency  of  the  coupled  system  are 
strongly  affected  by  the  vertical  advection,  and  the  larger  the  w  _li,  the  bigger  the  value  of 
£rr,  and  the  higher  the  frequency,  |  cr.  | . 

IX.  CONCLUSIONS 

A  coupled  cloud-OPBL  thermodynamic  model  has  been  presented.  This  model  demon¬ 
strates  the  positive  /  negative  feedback  mechanisms  between  clouds  and  ocean  mixed  layer, 
including  salinity  effects.  The  theory  predicts  the  generation  of  biennial  growing  oscillations 
at  typical  ranges  of  initial  states  in  the  tropical  ocean  and  atmosphere,  which  may  provide  a 
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new  explanation  for  the  low  frequency  oscillation  in  the  atmosphere  and  oceans. 

This  research  was  supported  by  the  Office  of  Naval  Research  and  the  U.S.  Naval  Postgraduate  School. 
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